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Abstract. Band structure trends in hole-doped cuprates and correlations with 
Tc are interpreted within the s-d exchange mechanism of high-Tc superconduc- 
tivity. The dependence of Tc on the position of the copper 4s level finds a natural 
explanation in the generic Cu 3d, Cu 4s, O 2pv and O 2py four-band model. The 
Cu 3d-Cu 4s intra-atomic exchange interaction is incorporated in the standard 
BCS scheme. This dependence of Tc in the whole interval of 25-125 K has no 
alternative explanation at present, and possibly this quarter of a century standing 
puzzle is already solved. 
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1 Introduction 

Quarter of a century after the dramatical discovery of high temperature su- 
perconductivity (HTS) the problem of its mechanism remains one of the 
longest-standing puzzles in the history of science. The cornucopia of ideas 
is immense - almost every quantum process was tested whether it is the 
long sought mechanism of HTS. In parallel, the intensive experimental re- 
search made cuprates into the most investigated materials. First-principles 
electronic structure calculations play an important role in understanding the 
physics of these materials. The Fermi surface of these highly anisotropic 
crystals is almost cylindrical with rounded-square cross- section. Fifteen years 
after the beginning of the cuprates era the eye of the professionalist has 
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uncovered a subtle correlation [1] between the shape of the Fermi contour 
and the critical temperature Tc for optimally hole doped cuprates. Pavarini 
et al. [1] considered many hole doped cuprates: Ca2Cu02Cl2, La2Cu04, 
Bi2Sr2Cu206, Tl2Ba2Cu06, Pb2Sr2Cu206, TlBaLaCuO,,, HgBa2Cu04, 
LaBa2Cu307, La2CaCu206, Pb2Sr2YCu308, YBa2Cu307, Tl2Ba2CaCu208, 
Tl2Ba2Ca2Cu30io, HgBa2Ca2Cu308, HgBa2CaCu206; the list can be ex- 
tended but we believe that a universal correlation was discovered. 

In an acceptable approximation the electronic band structure can be described 
by the four-band Linear Combination of Atomic Orbitals (LCAO) model with 
on-site energies Cs, Cp, ed> and hopping parameters igp, tpd, tpp- In this approx- 
imation the Hilbert space spans over the Cu 2>d^2_y2, Cu 4s, O Ip^^ and O 2py 
valence orbitals. The range parameter r = 1/2(1 + s) which determines the 
shape of the Fermi contour is determined by the LCAO parameters and the Fermi 
energy Ey 

s{Ef) = (e, - EF)iEF - ep)/{2t,pf. (1) 

The correlation between and r from the work by Pavarini et al. [1] is repro- 
duced in Fig. I. How such a correlation can possibly exist? Superconductivity 
is, certainly, created by some interaction, whilst the electronic band structure de- 
scribes the properties of independent electrons. Is it then possible to uncover an 
unknown mechanism of interaction investigating only the properties of nonin- 
teracting particles? The Tc-r correlation covers the whole temperature range of 
HTS, but now ten years it remains unexplained. We suppose that this correlation 
is as important for HTS, as was the isotope effect for the conventional phonon 
superconductors half a century ago. The aim of the present work is to interpret 
the correlation reported by Pavarini et al. [1] in the framework of some of the 
models of HTS. 



2 Tc-r correlation within tlie s-d theory 



Let us introduce a small modification of the s parameter, s = —^s. Its reciprocal 

tpd 

value 

s(e) (es -e)(e-ep)' 
has an energy denominator typical for the perturbation theory as applied to the 
secular equation of the generic four-band model 
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where Cp is the electron energy for the conducting d-band, Sx ~ 2sm{px/2), 
Sy = 2sin(pjy/2), and the wave functionis normalized ZJp + ^p+Xp + YjJ = 1. 
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Figure 1. Critical temperature of hole doped cuprates versus band structure parameter r 
by Pavarini et al. [1]. The theoretical curve is calculated according to Eq. (7) (see text for 
details). 



The approximate solution for small hopping amplitudes 
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has a simple interpretation: in the initial approximation we have a pure Cu 3d 
state with Dp « 1, in first approximation we have a linear dependence from 
the ipd amplitude on O 2p levels Xp, Yp oc ipd, and finally in second ap- 
proximation, the amplitude on the Cu 4s orbital Sp oc ipd^sp is included by the 
second virtual transition proportional to tgp. In short, the Cu 4s amplitude of the 
conduction band can be phrased as 3d-to^s-by-2p. 

As it was concluded by Pavarini et al. [1] that materials with lower eg tend to be 
those with higher observed values of Tc max- In the materials with higher Tc max 
the axial orbital is almost Cu 4s. What is the simplest interpretation of these 
observations? There is an emerging consensus that superconductivity in cuprates 
is created by some exchange interaction. The Pavarini et al [I] correlation gives 
that higher Tc max is determined by the highest Cu 4s amplitude Sp. Now we 
have to recall that the most usual 3d^s intra-atomic exchange has one of the 
largest amplitudes in condensed matter physics. We also have to point out that 
transition metal ion-ion exchange integrals Ju , which create antiferromagnetism 
in insulator phase are smaller than intra atomic exchange integral Jsd- 
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The intensive investigations of the Kondo effect have demonstrated that the anti- 
ferromagnetic sign of the two electrons s-d exchange is the rule and the ferro- 
magnetic sign is an exception. Incorporated in the BCS scheme for the calcula- 
tion of the order parameter S, the anti-ferromagnetic sign gives pairing in singlet 
channel with momentum dependent gap 

Ap(r) = s(r)xp, xp = SpDp. (5) 

This gap is included in the fermion excitation energy 

Ep = {'nl + Aiy/^, ^p=ep-£F, (6) 

and for the temperature dependent order parameter S(T) we have the standard 
BCS equation 

where (...) denotes momentum-space averaging over the Brillouin zone. For 
a pedagogical derivation of the BSC gap equation in the present notations see 
the textbook [2] and references therein. SUghtly below the critical temperature 
where the order parameter disappears S(Tc — 0) = the gap equation reads 



(8) 



Supposing that Jgd, being an intra-atomic process, is weakly temperature de- 
pendent we can determent its value using band parameters eg, Cd, Cp, isp, ipd, 
tpp, Ep and known Tc. Then we can calculate, for example, the dependence of 
Tc max by position of the Cu 4s levels eg. One illustrative example is shown in 
Fig. 2. The almost linear dependence has to have a simple qualitative interpreta- 
tion. Let us try to reveal this simplicity using the BCS interpolation formula 

fcBTc = 1.14?k^De-i/^(°)^. (9) 

In the present case wd is an energy parameter of order of the bandwidth, N{E-p) 
is electronic density of states per spin at the Fermi level and 

The interpolation BCS formula gives 

ln(r,) « ln(1.14wD) - s{Ep)/[N{Ep)J,a], (11) 

which according to Eq. (2) correlations reported by Pavarini et al. [1] are actu- 
ally correlations between the critical temperature Tc max and the BCS coupling 
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Figure 2. Logarithm of Tc (K) versus ecu4s. Science starts with simphcity - the almost 
linear behavior corresponds to the well-known BCS fomulae given by Eqs. (9) and (11). 

constant Jsd/s(i?F)- This general correlation is typically perturbed by stripes, 
inhomogeneities, magnetic phenomena, structural phase transition, apex oxy- 
gen, chemical substitution, and many other accessories of HTS cuprates. Nev- 
ertheless, they can be clearly seen for all hole doped cuprates. This qualitative 
agreement gives hope that four-band model with incorporated s-d exchange can 
become a standard model for superconductivity of cuprates. The key role of the 
s{Ef) parameter reveals why cuprates are unique for reaching high-Tc. Imagine 
that we can easily tune the position of the Cu 4s level eg. The pre-exponential 
factor wd is so big that we can easily reach a sauna-temperature superconductiv- 
ity if Eg is small enough. Maximal Tc is reached upon 3d-2p^s hybridization. 
For the Cu-O duet we have maximal triple coincidence of levels which ensures 
the success of the Cu02 plane. For all other combinations of transition metal 
with a chalcogenide s is much bigger. For 25 years many ways to decrease 
were empirically found: apex oxygen, bilayer hopping, pressure, etc. Perhaps 
only metastable artificial layers are not completely investigated. 

3 Computational method 

The parent Cu02 layer is an insulator with a half-filled conduction band. Doping 
with p holes per Cu atom results in metalization with hole filling f ~ \+P- The 
optimal doping corresponds to Pmax = 0.16 and /max = 0.66. We will use 66% 
hole filling for all examples in the present work. The Fermi level is determined 
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by the condition 

f={e{e^>Ep)), where 0(ep > i^p) = |^ tll<Z' ^^^^ 

The parameters of the four-band model can be determined by comparison with 
first-principles electronic structure calculations. For example, using the F point 
one can determine the on-site energies Cd = f-{px = 0,py = 0) for the conduc- 
tion band. In the present paper we use as an illustration a set of parameters (in 
eV) similar to Ref. [1] 

e, = 5.4, ep = -l, = 0, (13) 
^sp ~ 2, = 1.5, tpp = 0.2. 

eV Assuming max = 90 K for Cg = 5.4 we obtain, according to Eq. (8), 
Jsd = '2-4:4: eV. Then for so fixed Jsd we can calculate the correlations of Tc, max 
with the range parameter 

r = -7^ — r. (14) 
2(1 + s) 

Let us discuss briefly our finding. Within the four-band model one can derive an 
exact equation for constant energy contours 

-2f(e)[cos(p^) -I- cos(py)] -I- At'{e) cos(p^) cos(py) = ^(e), (15) 

where t{e), t'{e) and q{e) are polynomial functions of energy. The ratio of the 
effective intra-layer hopping parameters 

describes small variations of the shape of the Fermi contour influenced by the 
position of the Cu 4s level Cg. 



4 Conclusions 



Electronic structure experts should be proud that after many years of systematic 
research band calculations have revealed that Tc depends on the Cu 4s level. This 
numerical experiment is actually the crucial one for understanding the mecha- 
nism of HTS. The r parameter is introduced in electronic structure calculations 
in such a way, that for the first time we have a mechanism of a physical phe- 
nomenon possibly revealed by computer. 

We advocate a conventional theoretical explanation of the Tc-r correlations [1] 
which has no alternative among other theories of HTS. We have to wait for 
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the appearance of some other descriptions, because this important correlation 
covers a hole range of HTS. The LCAO model was only a tool for the theoretical 
analysis of the mechanism of HTS. It should be interesting to determine Sp 
and Dp directly from partial wave analysis at the muffin spheres directly by 
the electronic band calculations. The interpolating Xp = SpDp function can 
be directly substituted in the BCS gap equation (7) and the equation for Tc, 
Eq. (8). For the cuprates Xp well describes the experimental data for the gap 
anisotropy Ap on the Fermi surface. It remains to be seen if this mechanism is 
applicable only to cuprates. 

The iron-based pnictides are suspected to pursue an Oscar for supporting role. It 
will be extremely interesting to probe if Eq. (8) can describe the common trends 
for the Tc and gap anisotropy of ferro-pnictides. In this regards, finally we would 
like to re-analyze the "boost" that the physics of superconductivity got from the 
iron pnictides. Very recently, Mazin [3] recalled the Matthias' rules, well known 
in the physicists folklore: 1) A high symmetry is good; cubic symmetry is the 
best. 2) A high density of electronic states is good. 3) Stay away from oxygen. 
4) Stay away from magnetism. 5) Stay away from insulators. 6) Stay away from 
theorists. 

In order to emphasize the common properties of cuprates and iron pnictides we 
feel it compelling to rephrase these rules: (1) A high symmetry is good, square 
symmetry is the best, layered structure ensures empty s-band. (2) Having empty 
s-band, the Fermi level falls in the narrow d-band which ensures high density of 
states. (3) The oxygen and pnictide p-orbitals are perfect "go-between" for the 
transition metal 3d and 4s orbitals. These oxidants create significant 4s polariza- 
tion of the conduction 3d band. (4) We are just in the epicenter of magnetism; 
the s-d exchange amplitude which creates the ferromagnetism of iron, possibly 
creates the superconductivity of iron pnictides and cuprates. (5) Having narrow 
3d band in the case of half-filling we observe metal-insulator transition. Transi- 
tion to insulator phase is a hint for high density of states at the Fermi level for a 
doped compound. (6) Favorites of great socialists like Lenard, Kapica and Stark, 
have given similar advice. In order to surmount this weakness it is necessary to 
remedy the inferiority complex of the band calculators. One Nobel prize is a 
good initial dose of this treatment. Up to now there is no Nobel prize awarded 
to computational physics, but the considered in this paper T^-r correlations are 
a typical physical phenomenon whose nature is revealed by numerical calcula- 
tions. Prediction of an artificial structure with high s parameter, further synthesis 
by solid-state chemists and measurement of Tc may successfully initiate a new 
direction of intensive research. 

In conclusion, we briefly re-state the main property of high- Tc 
superconductivity — the triple coincidence of the d-, s-, and p-levels which 
ensures big enough BCS coupling constant. This triple coincidence is analogous 
to the sequence of vernal equinox, full moon, and Sunday (the end of the winter 
season, the end of the moon month and the end of the week): a triple holiday 
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creating the Great-day. In this sense the high-Tc superconductivity is the Easter 
of condensed matter physics. 

The authors are thankful to Prof. Ivan Zhelyazkov for a critical reading of the 
manuscript and to Alvaro de Rujula for his interest and comments during the 
scientific conference in memory of Matey Mateev, held in Sofia, 11-12 April 
2011. 
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